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Abstract. We study the effects of single impurities on the transmission in microwave 
realizations of the photonic Kronig-Penney model, consisting of arrays of Teflon pieces 
alternating with air spacings in a microwave guide. As only the first propagating mode is 
considered, the system is essentially one dimensional obeying the Helmholtz equation. We 
derive analytical closed form expressions from which the band structure, frequency of defect 
modes, and band profiles can be determined. These agree very well with experimental data 
for all types of single defects considered (e. g. interstitial, substitutional) and shows that our 
experimental set-up serves to explore some of the phenomena occurring in more sophisticated 
experiments. Conversely, based on the understanding provided by our formulas, information 
about the unknown impurity can be determined by simply observing certain features in the 
experimental data for the transmission. Further, our results are directly applicable to the 
closely related quantum ID Kronig-Penney model. 
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1. Introduction 

Research on photonic crystals, theoretical and experimental, has been sustained at a high 
intensity for several years, ever since the publications of Yablonovitch [1] and John 2J in 
1987. One interest for this research is the potentially very high number of applications in 
optoelectronics (see, e.g., Chapter VII of ji,). A substantial part of these studies concerns 
the study of impurities or defects in one, two, or three dimensional (ID, 2D or 3D) photonic 
crystals. 

Impurities or defects in photonic crystals may sometimes be unwanted but may also be 
extremely useful. For example, impurity states lying in a complete photonic band gap can be 
used for a waveguide and thus be an essential part of optical devices [4| . By introducing defects 
periodically in a perfect photonic array, coupled cavity waveguides are formed. The coupling 
of the cavity modes create impurity bands which have potential applications for the design of 
high-efficiency waveguides and waveguide bends [51 E] ■ It is important to mention that defects 
can be introduced in a controlled manner in photonic array experiments, see, e.g. [7]. There 
are promising theoretical results [8] that point defects, in particular a substitutional defect in 
3D crystals formed by a lattice of air spheres on a silicon background, may be used as micro- 
cavity for localizing light at a given point. Whence we appreciate the technological, as well as 
academic, importance of understanding the effects of various types of defects or impurities. 

There are several calculational methods for the treatment of impurities in photonic 
crystals. These are based on plane wave expansion of the fields ^] finite-difference time-domain 
algorithms [101 [2] ; variational methods [8] ; R-matrix methods [11] ; transfer-matrix methods 
|12j . combined, if necessary, with finite element methods; super-cell [13', Green-function, and 
tight-binding methods (see, e.g. Ref. [5]). Eigenfrequencies and eigcnfunctions of defects can 
also be calculated via the method of Sakoda et al. 14J based on the numerical simulation of the 
excitation process of the defect mode by a virtual oscillation dipole moment, in conjunction 
with the finite-time domain algorithm. 

On the other hand, for electronic periodic systems the calculation of impurity states dates 
back to the mid 50's with the introduction of effective mass theory [El [16]. Useful modern 
methods for the calculation of impurity levels are the super-cell tight-binding methods [T7] . 
which are applicable to shallow, deep and intermediate impurity levels [18j . 

In this work we investigate experimentally and theoretically the effect of single defects or 
impurities in the transmission of the electromagnetic field through arrays formed by Teflon 
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Figure 1. Experimental set-up. Left: Overview of the wave guide. Right: Enlargement of 
the part close to the electric dipole antenna and the two carbon absorbers. The top brass 
plate covering the wave guide is not shown. The wave guide has a total diameter of 78 cm, a 
depth A = 1 cm and a width B = 2 cm. The frequency range of the first propagating mode 
is from 7.5 to 15 GHz. The length of the Teflon pieces shown is d' = 4 cm and their index of 
refraction is n = \/2.08. In this picture, the spacing between all Teflon pieces is d = 4 cm. 

pieces alternating with air spacings. A closely related system to ours has been studied by 
Pradhan et al. [11] who looked for effects of isolated impurities in a system formed by an array 
of coaxial connectors. 

It is well known that point defects can produce localized states in the gaps [21 [Mj. The first 
experimental observation in photonics was made by McCall et. al. [20 . Can we observe these 
in our experiment? How else do impurities manifest themselves in the transmission curves? 
In this paper we investigate these questions for various types of single defects. We compare 
experimental results with those obtained by the transfer matrix calculations, equation (|13p . 
and point out the most prominent and typical transport features. 

One purpose of this paper is to show that our experimental set up can be used as a test- 
bed to study some of the phenomena occurring in more sophisticated (and expensive) photonic 
or electronic experimental arrangements. The second purpose is to show that our analytical 
expressions, derived in Section IV for the photonic Kronig-Penney model with single defects, 
are very helpful in elucidating the various features in the band structure; in particular, the 
frequencies of the defect modes and band profiles observed in the experiments. 

We remark that since our system is described by the one dimensional Helmholtz equation 
our results are directly applicable also to one-dimensional semiconductor (electronic) crystals 
formed, e.g., by sequences of quantum dots. Effects of irregularities such as an additional 
scatterer or a displaced quantum dot from its regular position have been the focus of many 
investigations, especially since these were observed experimentally by Kouwenhoven et al. [5T] 
in the case of electronic transport in heterostructures and by McCall [20] in the case of photonic 
crystals. 

2. Experimental set-up and the model 

The experimental set-up is shown in figure [TJ It consists of a brass ring where a microwave 
guide is cut out, in which several Teflon pieces, two carbon absorbers, and an antenna are 
inserted. Another antenna is fixed on the top plate (not shown in the figure) covering the 
waveguide. The antennas are connected to a network analyzer that allows the measurement of 
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Figure 2. The photonic Kronig-Penney model. The gray blocks are the idealized Teflon 
pieces of length d' , the air spacing is d. a„ and 6„ are the amplitudes of the wave function 
in the air spacing to the right of the nth Teflon piece. 



all matrix elements of the scattering matrix. The absorbers are used to eliminate reflection at 
the ends of the array and hence mimic a ID scattering system with open ends. This circular 
wave guide has been used, with metallic screws instead of Teflon pieces, to study the microwave 
realization of the Hofstadter butterfly and transport properties of ID on-site correlated 
disorder potentials [53]. 

In our experiment, the cut-off frequency for the lowest mode is /min = c/2B w 7.5 GHz 
and the second mode opens at 15 GHz. All results presented in this paper are restricted to the 
regime of the first propagating mode. Thus the system is effectively one-dimensional. Since 
the perimeter of the ring (« 234 cm) is much larger than the maximum wavelength used in 
our experiments (2 < A < 4 cm), our theoretical model assumes, as an approximation, a linear 
set-up (see figure [2]). 

For the lowest TE mode {E^ = 0) at / < 15 GHz the y component of the electric field Ey 
is also zero and the x component is 

E,{y, z, t) = Eo sin (^) $(z) exp(i2^/t). 

The wave function $(2;) obeys the Helmholtz equation 

where 



K m = 0, (1) 



kr = \r—k^nl--^ (2) 



/(27r/)2^, ^ 

and Ur = -/Jli is the position-dependent index of refraction. In the case the Teflon pieces are 
periodically spaced our system is the electromagnetic counterpart of the quantum ID Kronig- 
Penney model. 

We use the transfer matrix approach to calculate the transmission. Since both antennas 
are placed in the air, our transfer matrix should connect the wave functions amplitudes from 
air to air. Referring to figure [H the transfer matrix Q for a single cell connects the amplitudes 
(a„, 6„) and ( 
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and k (fc') denotes the wave vector in air (Teflon pieces) (c. f. equation Q). The transfer 
matrix for a single cell are 

Qii = [cos(fc'd') + ia+ sm{k'd')] e''"^ (4) 

(5i2 = — za_ sin(fc'(i'), (5) 

where 

- '«> 

Due to preservation of flow and time reversibility it yields Qn = Q22, Q12 = Q21 and 
det{Q) = 1. The elements Qij are the same as in the quantum ID model of a square potential 
well, but with k = y/%lE, k' = ^j2^i{E - Vb) and Vb > 0. From (2) k' > k; thus, the Teflon 
piece acts as a well (Vb < 0) in quantum mechanics, except that in the photonic array the 
"depth" increases with frequency. 

The microwave vector network analyzers measure the full scattering matrix S defined by 

In terms of the Q matrix elements, the S matrix reads 

_Q21 

1 

The transmission T through a single cell is given by 

T=|5i2p = ^ = -4^> (9) 




(8) 



IQ22P 1 + IQ 



12 



with Q12 and Q22 given by Eqs. (4) and (5) for a single cell. 

Expression (9) together with equation ^ shows immediately that a single Teflon piece 
becomes completely transparent at frequencies obeying the relation 

k'd' = rmT, m= 1,2,3,... (10) 

We shall refer to these frequencies as 1-Teflon resonances. 

For an array of N equally spaced cells without any impurities or defects, we need to 
evaluate Q to the Nth power; a numerical process that can be easily performed. However, it 
is more illustrative to use the Cayley-Hamilton theorem of linear algebra to get [Ml HSl 



where 9 is the Bloch phase corresponding to the infinitely periodic array: 

cos 6' — "SiQii — cos{k' d') cos{kd) — a+sm{k'd')sm{kd). (12) 

Note that T = 1 not only when Q12 = (i. e., at the 1-Teflon resonances) but also whenever 
N9 = ztmr, n = 1, 2, . . . , iV — 1. The latter one gives rise to — 1 peaks with T = 1 in each 
transmission band since 9 shifts through tt [24l |25l [26l EI] ■ According to the above condition, 
the system of A'^ = 2 Teflon pieces separated by air becomes transparent when 9 — 7r/2, which 
in turn implies TrQ = 0. We call these the 2-cells resonances emphasizing that it is not just 
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two Teflon pieces next to each other but separated by an air spacing. A httle thought reveals 
that the — 1 oscillations are centered around the 2-cells resonance frequency (see also [25]). 

Luna-Acosta et al [53] treat the case of the regular Kronig-Penney model (i.e., periodic 
and no defects) where it is shown that equation (jlip reproduces most details of the experimental 
transmission data as a function of frequency for all kinds of periodic arrays considered. Since 
we are concerned here with the effects of the impurities, we show for comparison the theoretical 
and experimental results for the transmission amplitude |>S'i2| for an array 16 equally spaced 
cells, with Teflon pieces of width (i'=4.0cm and air spacing of width d= 4.0 cm, see flgurelDJa). 

For the experimental and theoretical curves to agree as well, it was necessary to define 
some effective length for the length of the Teflon pieces and the air spacings. That is, since 
the actual waveguide is circular, the Teflon pieces and air spacings are not rectangular pieces 
but slightly curved with the larger side being exactly 5 percent larger than the shorter side. 
(The 4 cm Teflon pieces and the air spacings referred to above actually mean that the shorter 
side is 4 cm whereas the larger one is 4.2 cm.) It turned out that the best fit could be obtained 
with an effective length of 4.08 cm, i.e, a 2 percent larger than its shorter side. We emphasize 
that this is the only fitting parameter in our calculations and moreover it is the same for all 
calculations presented here. Throughout the paper, whenever we quote the width of a Teflon 
piece or an air spacing we mean the shorter side of it and in the numerical calculations we use 
their corresponding effective length. 

We remark that our model does not consider the absorption of the signal. However, 
comparison of the experimental and theoretical data shows that the band structure (gaps and 
band profiles) is not affected by the absorption except for the attenuation of the transmission, 
which is about constant throughout the frequency range (the experimental transmission is 
about five times weaker). Thus, absorption in our experiment, does not destroy coherent 
phenomena hke the band structure (see also [23 )• 

In figure [3] and all subsequent transmission plots the frequency values of the 1-Teflon and 
2-cells resonances are marked by crosses and circles, respectively. Different types of bands are 
formed depending on the position of the Teflon resonances relative to the 2-cells resonances 
[25j . In these and all forthcoming transmission plots we mark with shaded strips the gap 
regions, defined by the condition that the eigenvalues X± of the transfer matrix Q are real. 
Since A± = <^ ± a/^^ — 1, where ^ — ^Qii, forbidden gaps occur whenever > 1. 

Given the good agreement between experimental transmission data and the photonic 
Kronig-Penney model, implied by equation (llip , we proceed to discuss the effects of impurities 
in the photonic Kronig-Penney model. 

3. Impurity eflFects in the transmission 

Different types of defects or impurities can be realized in our experimental set-up with just 
air and Teflon segments and can be represented by one of the two general sketches shown in 
figure m The upper sketch illustrates a general interstitial impurity: an extra piece of Teflon 
of width b placed somewhere in the air spacing between two pieces of Teflon. The rest of the 
arrangement remains unaltered; the perpendicular dotted lines mark the boundaries of the unit 
cell along the pure crystal. By definition of point defect, the interstitial impurity breaks the 
periodicity only locally. Hence, the length of the extra piece of Teflon should be less than or 
equal to the air spacing. If the width b is larger than this, then the boundaries of the regular 
cells would be displaced from their original position and consequently periodicity would be 
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Figure 3. Interstitial impurities for a 16 cells array with d = d' = 4 cm. Dotted (solid) 
lines are the experimental (theoretical) curves. The shaded regions mark the forbidden gaps, 
given by the condition < 1 (a) pure (impurity- free) array; (b) small Teflon piece (3.16 cm) 
inserted in the center between the 3rd and 4th regular Teflon pieces, (c) small Teflon piece 
(3.16 cm) inserted in the center between the 8th and 9th regular Teflon pieces. The crosses 
(circles) mark the frequency values for the 1-Tefion-resonances (2-cells resonance), see text. 
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Figure 4. Sketch of impurity arrangements. In the upper part the interstitial, in the lower 
part a substitutional impurity is shown. In the case of the substitutional A^i + N2 = N — 3, 
whereas in the case of the interstitial A''i + A'2 = N — 2, where is the total number of 
unit cells of the perfect crystal. 



globally broken (it becomes an extended defect, e.g., a topological defect [551 [IS])- In this 
paper we do not consider topological defects. 

As shown in the upper sketch there are A^i unit cells periodically arranged to the left of 
point xi, each described by the matrix Q, then periodicity is interrupted by the impurity at 
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point xi and again recovered at point X2. To the right of point X2 there are N2 cells periodically 
arranged. The total transfer matrix going from left end of the array to right end of the arrays 
can be written as 

Qtot — Q^^QimpQ^^ , (13) 

where Qimp is the transfer matrix connecting the amplitudes of the wave function at xi with 
those at point X2- Specifically, 

Qimp = Dd/2MDcMtDaMDd/2, (14) 

where Di (i = d/2, a, and c), is the transfer matrix corresponding to an air spacing, of length 
i. M is the transfer matrix for the regular size Teflon piece and Mi, that of the Teflon piece of 
length 6. Note that Ni + N2 = N — 2, where N is the total number of unit cells (and Teflon 
pieces) of the perfect crystal. 

The general case of a single substitutional impurity is shown in the lower sketch of figure 3) 
The substitutional impurity has a length b, shown here longer than the regular size Teflon piece, 
but it could be of any length and/or of a different dielectric material. Quick inspection shows 
that the form of the impurity matrix is the same, but the values of the lengths a, b, and c are 
different and now A^i + A^2 = N — 3. 

Thus, we have a straight-forward numerical scheme to calculate the transmission amplitude 
in the presence of an impurity: for a given impurity arrangement we compute equation (|13p . 

In figure ^h) we show the case when an extra Teflon piece of smaller length (than the 
length of the default Teflon piece) is inserted in the center of the air spacing between the 3rd 
and 4th regular Teflon pieces, thus iVi = 2,N2 = 12). Figure [3jc) shows the case of the same 
type of impurity but located in the air spacing between the 8th and 9th regular Teflon pieces, 
thus Ni = N2 = 7. Notice that the good agreement between the transfer matrix calculations 
and main features of the experimental data. 

Inspection of the plots of figure [3] reveals that the impurity affects the transmission in 
two ways: i) The bands develop different profiles (compared to the pure array, figure[3ja)). ii) 
Peaks appear in most of the gaps in figure 3(b) coinciding in frequency with those of figure[3ljc). 
Such a coincidence is expected since the impurity in both arrays occupies the same relative 
position within the cell (in the center of the air spacing) but in a different cell. 

In figure [3Jb) the band profiles look like a superposition of slow and fast oscillations. 
In figure [3Jc) the oscillations are larger and only of one type. In both, figure [3Jb) and (c), 
bands containing 1-Teflon resonances (see the 2nd, 4th, and 6th bands) are less affected by the 
impurity. Note also that the intensity of the peaks is greater when the impurity is placed near 
the center of the array, figure [3fc), than when it is placed near the end of the array, figure [3fb). 
This feature was invariably found in all our experiments. 

Figure [S] shows the experimental and theoretical curves for three types of substitutional 
impurities. Figure [5fa) pertains to the case when a Tefion is removed from the array, i. e., a 
vacancy. Specifically, the 8th Teflon is removed. In figure [DJb) the 8th regular size Teflon piece 
is substituted by a smaller, whereas in (c) it is substituted by a larger one. Note that while 
the defect is in the same number of cell within the array, the band profiles, the frequencies and 
intensities of the peaks are different for each case since the defects are different. 

A single point defect can produce multiple impurity levels within a band, as has been 
observed, e. g. by Yablonovitch et al. jSOj in an three dimensional photonic array with a donor- 
like impurity. Note that the figure [5l^c) shows also two spikes in most of the gaps. 



CONTENTS 



a) 



b) 



c) 



10° 
10-1 

10-2 
10-3 

10° 
10-1 

10-2 
10-3 

10° 
10-1 

10-2 
10-3 



i I 




V 








- 


- \ 






/Arv~^---^x\/Y 


1, ^ 


i 


i • 


"r i 














1 


1, ^ 









10 



1 1 



14 



15 



16 



12 13 
f [GHz] 

Figure 5. Three types of substitutional impurities in an 16 cell array {d = d' = 4cm). 
Dotted (solid) lines are the experimental (theoretical) curves. In (a) the 8th Teflon piece 
was removed corresponding to a vacancy. In (b), (c) it was replaced by a smaller (larger) 
size Teflon of length d' = 3.16cm {d' = 6.32cm), respectively. 



Shallow impurities in semiconductors, associated with long-range binding, typically show 
up as spikes near the band edges, whereas deep impurities typically lie near the center of the 
gap [18]. However, figure [5] shows that a single impurity may produce spikes near the center of 
the gap or near the band edges. Thus we see that the above characterization of the impurity 
as either deep or shallow is not always applicable. 

Although our numerical procedure is seen to reproduce quite well the main features of 
the experimental data, it does not serve to explain their origin. For this purpose we develop 
a scheme in the next section which provides us with a general understanding of the direct 
relation between the number of slow and large or fast and small oscillations and the position 
of the impurities in the array. We will also show how to determine the frequency position of 
the impurities. 

4. Analysis of the effects of the impurities. 

We express the transfer matrix Q in terms of its diagonal representation A. Inserting 
Q PAP^^ into equation ([14|) to get 

Qtot = PA^^gA^ip-i, (15) 

Q = P-^Qin.pP, (16) 



where 



with P = { ^--Q" ) , A = ( ^+ ^ ) , A± = e ± TF^, and ^ = TrQ/2. 
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Note that Q is the impurity matrix in the representation of the diagonal basis of Q. This 
"rotated impurity matrix" contains information about the coupling of the impurity to the host 
material through the elements of the transformation matrix P. 

From the definition of Qtot and using the representation of P given above we get after 
some algebra useful expressions for the matrix elements of Qtot^ 

Qtot,22 = A^'^+"^ (/+0n + Z+g^iA;^^^ - AQi^A;^^^ - Ag22A;'(^^+^^)) , (17) 
where f± = (A± - Qii)/{X+ - A_), 

Q21 = Qi„.p,ii/- + %Mi(A_ - Qn)/+ - - Qi„,p.22/-, (18) 

V12 ^+ — Wii 

and 

<922 = Qimp,ll/+ + ^'"'P-^^ _ Qii)/+ — ^21IlEl^lRll — Qiinp 22/-- (19) 

tt/12 •A+ — yii 

The terms Qn and Q12 are obtained from Q22 and Q21, respectively, by exchanging A+ with 
A_ everywhere. As in the case of a single scatterer, it is also true that the complex conjugate 
of Qtot.ii ((3tot,i2) is (9tot,22 ((3tot,2i)- This property is not obeyed by the elements of the 
rotated impurity matrix. On the other hand, because equation (|16p is a unitary transformation, 
Tr(Q) = Tr(gi^p) and det(g) = det(gi^p). 

Setting Ni + N2 = iV, Qn = Q22, and Q12 = Q21 = in equation ((T6l) and using 
T = l/|(5tot.22| we recover equation (jlip that is valid only for the pure array. 

4-0.1. Band oscillations To analyze the effects of the impurity on the band profiles it is 
convenient, since T — |Qtot,22|^^j to study in detail the element Qtot.22 given by equation (|17p . 
In the bands the eigenvalues are complex, A_ = A!^ and of modulus one, so we can (when 
considering the absolute value of (5tot,22) ignore from our analysis the factor x^^+^^ front of 
expression (17). We now discuss the quantities f± and Qij. The functions f± do not contain 
any information on the impurity: they depend only on A± and Qii, characterizing the periodic 
array. f± are real and smooth ( almost flat) functions within the bands, diverging as l/-y/l— ^ 
at the band edges [J. More important for our analysis is to note that /+ and /_ take alternate 
roles in consecutive bands. Namely, |/+| is larger than |/_| in the first and all odd- numbered 
bands, while the opposite is true in the second and all even-numbered bands. 

The quantities Qij are also typically very smooth functions within the bands, diverging 
like — ^"^ at the band edges, and are all roughly of the same magnitude, except for Q21 

which diverges at the Teflon resonance (but note that /+ = at Teflon resonances). 

With this information about the f± and Qij it is expected that in the odd-numbered bands, 
where terms containing /_ can be ignored to first order, the first two terms in equation (jl7p give 
the main contribution to Qtot,22- Specifically, T « |/+(Qii + Q2iX+^'^^)\'^ should yield N2 
oscillations (recall that 9 shifts through tt in each band). Figure [SU a) shows this term and also 
as a function of frequency for the case of an interstitial impurity centered between 

f In the bands: \± = g it 17 with 7 = — real, and Qn can be written as Qii = i; + irj. It follows 
from the definition of f± that /„ = — {7 + ?7)/27 and = /_ + !. Given that 77 and 7 are smooth functions 
everywhere and 7 goes to zero at the band edges (where = 1) then f± are very smooth functions within the 
bands, diverging as ~ at the band edges. 
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Figure 6. Band oscillations due to impurities, (a) and (b) correspond to an interstitial 
impurity centered in the air spacing between the third and the fourth regular size Teflon 
pieces, whereas in (c) it is displaced. |/_|_Qii|^^ (dashed) and |/_|.Qii + f+Q2i^^'^^^\~^ 
(solid) are shown in (a). In (b) and (c) |Qtot,22|^^ (dashed) and |/-|_Qii — f—Qi2^^^^^ 
(solid) are shown. 



the third and fourth Teflon piece (the case corresponding to figure [31(b)). For the moment 
let us concentrate on the first and third transmission band. The term |/+Qii|^^ (dashed 
curve in (a)) gives the zero order approximation to the transmission band. According to the 
notation of figure |4] it foUows that Ni = 2 and N2 = 12. The 12 oscillations can be counted in 
figure SKa) in the plot of |/+(Qii + Q2iXf^')\-^. Moreover, since |/+(Qii + Q2iX+^^')\-^ 
gives the dominant contribution in the odd-numbered bands, the 12 oscillations should also be 
observed in the total transmission. Indeed this is the case, as it is shown by the dashed curve 
in figure [Sfb). This plot also shows that the total transmission curve oscillates about a slowly 
oscillating curve (thick solid line) with two maxima. This slowly oscillating curve is in fact 
given by \ f+Qii — /-Qi2AI^^^ with iVi = 2 as pointed out above. Finally the last term in 
equation (fT7|) gives additional A^i + N2 oscillations of smaller amplitude, which provide small 
corrections (not readily visible in the plot). 

As figure [6jb) shows, the term \ f+Qii — f-Qi2^Z^^^\^^ does not agree at all with the 
average curve of the transmission in the gaps nor in the even numbered bands. This is 
because our discussion has been limited to the bands, where the eigenvalues are complex, 
and of modulus one, whereas in the gaps the eigenvalues are real, and hence the factor 
^Ni+N2 J- 2 7^ cannot be ignored (in the next section we shall consider in detail the gap 
regions). There is no agreement in the even- numbered bands either because there the dominant 
terms are those containing the form factor namely, the sum |/_((522 + Qi2'^I^^^)|- 
This can be readily seen by factoring in equation (fTT]) the term ^^^1+^2 Qtot,22 = 

-A^^+^^(/_Q22 + /-Qi2A:2^^ -Z+QsiAl^^^ -/+QiiA:^^^'+^"'). written in this form, it is 
clear that the same arguments used for the odd-numbered bands work for the even numbered 
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bands but with /_ and interchanged. Exphcitly, the dominant term in the even-numbered 
bands is |/_(Q22 + 02i^I^^^)|7 giving rise again to N2 — 12 oscillations. The second order 
term, producing the Ni — 2 oscillations, is I/-Q22 — /+Q2iAl^^^|. 

Thus, the number of slow and fast oscillations appearing in the profiles of the transmission 
bands can be read off directly from equation |J7[ j. Assuming for the moment that N2 > Ni, the 
band profiles should show N2 fast oscillations about a curve with iVi slow oscillations. This is 
exactly what we have noticed in our experimental and numerical transmission curves, see, e. g., 
figure [Sljb), corresponding to an interstitial impurity with {Ni,N2) = (2, 12), respectively. 

On the other hand, if Ni = N2, then the profile should show only Ni oscillations of larger 
amplitude than in the case above (iV2 ^ Ni). The experiments confirm this, as it is illustrated 
in figure [2^c), where Ni — N2 — 7. Finally, if A^i is close but not exactly equal to iV2 then 
the oscillation pattern becomes more complicated, with irregular (incomplete) oscillations and 
with no distinctive modulation pattern, as exemplified by the plots of figure [5l where Ni=6 
and iV2=8. 

Thus, given some experimental curves for the transmission, one can determine where along 
the array there is a single impurity by simply observing the pattern and number of oscillations. 
Note, however, that due to time reversal symmetry it is not possible to deduce on which side of 
the array to count the A^i cells. At the very least, the form of the profile indicates the presence 
of impurity or defect, and whether it is located near the ends or near the center of the array. 

It is important to remark that the above conclusions hold as long as the elements Qij do 
not have oscillations of their own as it is the case of the interstitial impurities discussed above. 
For certain other types of defects the elements Qij do have some structure and consequently 
the type of profiles discussed above become somewhat distorted. An extreme example occurs 
when a smaller Teflon piece is placed off-center between two regular size Tefion pieces, that is, 
a displaced interstitial impurity. Specifically, Qij show a strong frequency dependence (e.g., 
Qii and Q21 have each three extrema in their real and imaginary parts in each band). For 
example, figure[6l^c) shows the total transmission and jJ+Qn — f-Qi2^^^ for an interstitial 
impurity placed off-center between the 3rd and 4th regular size Teflon. This defect is described 
by the upper sketch of figure [4] with A^i = 2 and N2 = 12 with N = 16. Note that the term 
IZ+Qii — f-Qi2^-^^ correctly gives the average curve of the band profile but now there are 
three slow oscillations and 14 fast oscillations, instead of 2 and 12 oscillations, respectively. 

4-0.2. Impurity states We now analyze the appearance of the peaks in the gaps, the well 
known signal for the presence of an impurity. Without limiting in any way the conclusions to 
be drawn here, let us consider the d = d' =4 cm array, the example we have been using, and 
let us label the gaps. Gap one is in the frequency range [7.7, 8.1] GHz, see, e.g., figure [TJ gap 
two in the range [8.6,8.9] GHz, and so on. 

Recall that in the gaps the eigenvalues are real, therefore, we rewrite equation (|17p as 

Qtot.22 = /+QiiA^ + /+02iA^^"^^ - /-Qi2A^^-^^ - /-Q22A^, (20) 

where N — Ni + N2. Now, for the transmission to produce a sharp peak in the gaps, 
Qtot,22 must be small in a very narrow range of frequencies. Note that A+A_ = 1, with 
A_ and A+ alternatingly being the larger one in each gap. In the even-numbered gaps (where 
A-|_ > 1 > A_) the dominant term in (5tot.22 is the first term in equation (|20p . and it will be 
large (corresponding to a negligible transmission) unless it happens that Qn is very small or 
zero. Similarly, in the odd-numbered gaps (where A_ > 1 > A+), the dominant term is the 
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Figure 7. Impurity states. In both plots the transmission amplitude l/|Qtot,22| is shown 
(solid). In (a) the dashed lines corresponds to Q22 and in (b) to Qii. All parameters are 
the same as in figure ISfc), except now the impurity is between the 7th and 8th regular size 
Teflon piece, thus A''i = 6 and N2 = 8. 

last one and it will also be a large number unless Q22 happens to be very small. To illustrate 
this general observation, in figure [7| we consider the case of an interstitial impurity placed in 
the air spacing of the 8th cell in a 16-cell array (thus, A''! = 6,A^2 = 8). In figure [Tfa) we 
show the transmission amplitude l/|(5tot,22| and also the term |Q22|- In contrast, in (b) the 
term |Qii| is shown, together with l/|(5tot,22|- Observe in figure[7l^a) that whenever Q22 goes 
sharply to zero in the odd-numbered gaps (A_ > 1), there is an impurity state at the same 
frequency. Q22 also goes to zero sharply in the even-numbered gaps(A+ > 1) but here there 
is no coincidence with the impurity states, in agreement with our argument of the previous 
paragraph. Similarly, figure [71[b) shows that there is an impurity state in the even-numbered 
gaps at the frequency values where Qn goes to zero. Thus, our simple argument above really 
works in determining the position of the impurity states: an impurity state occurs in odd- 
numbered gaps at the frequencies where Q22 'is zero and in the even-numbered gaps where Qn 
is zero. A similar case was already shown in figure ^ic) except there (A^i = N2 = 7 and 
thus band profiles are different but peak positions are the same), where we see an excellent 
agreement with experiment concerning the location of the peaks. 

Single impurity states are known to decay exponentially away from the site of the impurity. 
One then may expect that the intensity of the peak should be stronger the closer the impurity 
is from any of the ends of the array. However close inspection of formula (20) indicates 
that this is not the case. That is, exactly at the impurity state frequency, say when Qn 
is zero in an even-numbered band, the first term in (20) is zero and the leading terms are 
f+Q2i^+^~^^ — /-Q12 . The inverse of its absolute value gives to a good approximation 

the intensity of the peak. For simplicity, let us consider first the case /+(52i = /-Q12 which 
is valid only for "symmetric impurities", i.e., when the distance a equals the distance c (see 
figure m). Note that ^JV^-JVi _^^JV2-Ari 

as a function of N2 — Ni is symmetric about its minimum 
N2 — Ni — 0. Thus, recalling that the transmission is the reciprocal of |Qtot,22|, we conclude 
that the intensity of the impurity state should be stronger the closer the impurity is to the 
center of the array. It is strongest when N2 ~ Ni, i. e., when the impurity or defect is at the 
center of the array. 
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The above argument assumed that the impurity is symmetric. If it is not the case, then 
it can be shown that the minimum of \f+Q2i^+^~^^ — f-Qi2^^^^^^\ does not occur when 
the impurity is at the center (A^i = A*'2) but at a distance A^i — N2 that is proportional to 
the logarithm of |/_(5i2|/|/+<92i|- Since this ratio is not large, the highest intensity peak 
occurs still when the impurity is close to the center of the array. This was verified in all our 
experiments and can be seen clearly by comparing figure [3Jb) and figure [3{c) . We emphasize 
that this effect is independent of the type of the periodic potential and type of impurity in 
a ID periodic array. In fact, the same behavior was noted in an experiment with a coaxial 
connector photonic crystal [19j where, however, no explanation was given. 

As we have seen, the intensity of the peaks always decreases as the number N of 
cells increases. However, an important conclusion, drawn from a detailed examination of 
equation (I20|) as a function of A+, for fixed N, is the following: the decrement in the peak 
intensity as N increases is weaker the farther the impurity state is from the center of the 
transmission band. This was confirmed by our experimental and theoretical results and it is 
in agreement with the generally accepted idea that shallow states near the band edges are the 
shallow levels, associated with long range potentials, in contrast with the deep levels lying near 
the center of the band (see, e.g. .[15]). 

5. Conclusions 

We analyzed the effects on the transmission of single impurities in an array of regularly spaced 
pieces of Teflon in a microwave guide, described by a ID photonic Kronig-Penney Model. We 
performed a series of experiments with point defects of various types; namely, interstitial, 
substitutional, and vacancies and showed that single impurity affects the transmission in two 
ways. One is the well known appearance of localized states in the gaps; the second, to our 
knowledge not discussed so far in the literature, is the appearance of fast and slow oscillations 
in the band profiles. Our transfer matrix calculations correctly predicted these features. 
Experiments with other types of single impurities, e.g., displaced interstitial and displaced 
substitutional impurities, not reported here, were also found to be equally well described by 
our model. The transfer matrix calculations (being purely numerical) agree very well with the 
experimental data, but they provide no insight for the understanding of the various features 
observed. Thus, as an important contribution, we derived an exact closed form expression 
for the transmission amplitude that is useful in elucidating the effects of a single impurity on 
the band profiles, the impurity levels, and their intensity. This expression involves elements 
of the transfer matrix Q of the regular cells in the array and of the transfer matrix Q of 
the impurity. Since this matrix is written in the representation that diagonalize the transfer 
matrix of the regular cells, Q contains information about the impurity and its coupling with 
the host environment. We found that the set of impurity levels is given by the zeroes of the 
diagonal elements of Q. Further it was shown that the intensity of the impurity states in the 
gaps depends on two factors, namely, the off-diagonal elements of Q and the distance of the 
impurity from the center of the array. The closer the impurity is to any of the ends of the array, 
the lower the intensity of its level, and vice versa. It was also shown that the number of fast 
and slow oscillations in the bands give direct information about the position of the impurity, 
relative to the center of the array. 

We note that the agreement between the experimental results and analytical calculations 
gives us confidence to treat the inverse problem. That is, we can extract information about 



CONTENTS 



15 



the unknown defect from the inspection of the band profiles and the locahzed defect modes. 

Our method is an alternative approach to the tight-binding and Green's functions methods, for 
ID systems, with the advantage that it is simpler and elucidating. Although the experimental 
realization here is in the microwave regime, the model and the results are equally valid for 
higher scales of frequency corresponding to light experiments. Wc remark that we can apply 
our formalism in a straight-forward manner to any kind of regular array with single defects once 
we have determined the particular transfer matrices for the regular cell and for the impurity. 

Wc emphasize that the formulas (17-20) and the procedure for determining the position 
and intensity of impurity states and features of the band profiles for the photonic Kronig- 
Penney model are exactly the same as for the electronic Kronig-Penney model (with square 
barriers) for energies above the barrier. This is so since their transfer matrix elements are 
identical and the difference is only in the definition of the longitudinal wave vector. Hence 
our procedure is useful for analyzing the effects of impurities either, in the transport of 
charged particles in ID periodic arrays of electronic potentials (e.g., heterostructures) or the 
transmission of electromagnetic waves in ID photonic arrays. 

Finally, the fact that our experimental and theoretical results agree very well give us 
confidence to treat various other types of specific arrangements with the goal of realizing them 
in other experimental set-ups as photonic devices. We would like to stress that this is by 
no means self evident. The observed transmission patterns are the result of a complicated 
interplay of interferences from all structures within the waveguide, and it is well known 
that interferences are extremely sensitive to perturbations, in particular in the presence of 
absorption. The present work has shown that we can rely on the experiment in this respect. 
Thus our experimental set-up can be used as a preliminary and low cost study of a 
particular Kronig-Penney model to be implemented later in more sophisticated and expensive 
microscopic photonic realizations, which might have real applications. 
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